Surface Region of Superfluid Helium as an Inhomogeneous Bose-Condensed Gas 
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We present arguments that the low density surface region of self-bounded superfluid 4 He systems is 
an inhomogeneous dilute Bose gas, with almost all of the atoms occupying the same single-particle 
state at T — 0. Numerical evidence for this complete Bose-Einstein condensation was first given 
by the many-body variational calculations of 4 He droplets by Lewart, Pandharipande and Pieper 
in 1988. We show that the low density surface region can be treated rigorously using a generalized 
Gross-Pitaevskii equation for the Bose order parameter. 
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PACS numbers: 05.30.Jp, 67.40.-w, 03.75Fi 

In recent years, there has been a major effort to pro- 
duce atomic-like gases in a Bose-condensed phase, using 
excitons, spin-polarized Hydrogen and laser-trapped al- 
kali atoms (for reviews of this work, see Refs. In 
such systems, the fraction of atoms in the lowest energy 
single-particle state could approach 100% at T = 0. This 
unique phase of matter was first studied by Einstein and 
London and would exhibit many unusual properties [Q, 
including superfluidity. In contrast, in bulk superfluid 
4 He, the condensate fraction is only about 10% at T = 
Thus the ground state wavefunction does not have 
the simplicity of a gas with complete Bose-Einstein con- 
densation (BEC). 

In the present letter, we argue that such a dilute Bose- 
condensed gas is in fact already present in the surface 
region of superfluid 4 He, where the density becomes very 
small. This striking phenomenon has apparently been 
overlooked in the extensive theoretical literature on sur- 
faces, films and droplets of superfluid 4 He (for references, 
see Ref. 0). Moreover, in this low density surface region, 
where one has almost 100% BEC into a single state, one 
can do microscopic calculations of the order parameter 
using the standard field-theoretic description of an inho- 
mogeneous Bose-condensed fluid [|[^| . Our discussion is 
limited to ground state properties (T — 0) but we make 
a few remarks about finite temperature at the end. 

Striking confirmation of the preceding physical argu- 
ment is available in the results of Lewart, Pandharipande 
and Pieper (LPP) M, who studied 4 He droplets (with up 
to 240 atoms) using a variational ground-state wavefunc- 
tion approach g. These authors calculated the single- 
particle "natural orbitals" ipi( r ) which diagonalize the 
single-particle density matrix, giving 



Pl (r,r')=]>>^*(r)^(r'), 



(1) 



where n, is the occupation probability of the state ipi( r )- 
The local density is given by p(r) = pi(r, r). By direct 
numerical calculation, LPP found that the node-less ls- 
state ipi s (r) had an occupation probability of n ls ~ 25, 
while for all the other states, rij < 0.5. They identified 
(and justified) this Is state as the Bose-condensed state, 
with an associated local condensate density being given 
by 



p c (r) = ni s |^i s (r)| 5 



(2) 



Their results for 
However, LPP did 



LPP evaluated the total local density p(r) and the con 
densate local density p c (r) using (2) 
N = 70 atoms are replotted in Fig. 1. 
not explicitly discuss the implications of the most striking 
feature shown in Fig. [l], namely that since in the surface 
region almost all the atoms at a given value of r are in 
the ipisi?) state, this surface region provides a physical 
realization of the long sought-for dilute Bose-condensed 
gas §!• 

An important aspect of the calculations in Ref. j7j is 
that near the center of the drop, p(r) and p c ( r ) are very 
close to the values in bulk superfluid 4 He, with p c (r) ~ 
0.1p(r). This indicates that even droplets of N = 70 
atoms are large enough to ensure that the surface region 
density profile will be very similar to that of a planar 
free surface of bulk superfluid 4 He. This is confirmed by 
more recent calculations on larger droplets ||. In the 
case of planar (vs spherical) symmetry, the atoms will 
Bose condense into states with zero momentum parallel 
to the surface [id . 

Krotscheck [M first found that the local condensate 
density p c ( r ) increased when the total density p(r) de- 
creased in liquid 4 He in contact with a wall. Camp- 
bell flll|] h as recently drawn attention to the results in 
Ref. [|7| -pi and discussed them in the context of the for- 
mal theory of off-diagonal long range order (ODLRO) 
in liquid He droplets. Apart from these papers, Bose- 
condensation has been ignored in the literature on inho- 
mogeneous Bose liquids at T — (for references, see Q] 
and p^]). Moreover, even in Refs. ]7|-pd]], what is still 
missing is a physical interpretation of the numerical re- 
sults for p c (r) at surfaces as evidence for (and relevance 
of) a low density Bose gas in which the condensate frac- 
tion can be 100%, and the important implications of this 
new picture. It is this aspect which the present letter 
addresses. 

The numerical results in Fig. |l| are consistent with the 
empirical formula given in Ref. 



Pc(r) = p(r) 



1-0.68 



p(r) 



Pb 



(3) 



where pb is the bulk or saturation density of superfluid 



1 



4 He at T = 0. It is easy to verify using (^) that the 
local condensate fraction p c (r) / 1 p(r) smoothly increases 
from 0.1 to unity as p(r) goes from ps to zero. One 
of the most interesting features shown by Fig. [I] is that 
p c (r) develops a maximum. Using ([}]), this occurs when 
p{r) ~ 0.5ps and corresponds to p c (r) — 0.22pB- Thus 
the results of Ref. JjJ predict that in a certain region, the 
density of 4 He atoms which are Bose-condensed can be 
twice as large as the bulk value p c ~ O.lps! Of course, as 
p(r) decreases to zero, so must p c {r). While the precise 
magnitude may depend on the approximate calculations 
used in 0, we believe that the increase in p c (r) as we 
go from the center of the 4 He drop is a real effect. As 
noted in Ref. |Q, the initial increase in p c (r) shown in 
Fig. |l| for p(r) < ps is consistent with the calculated and 
observed decrease in p c in bulk liquid 4 He under pressure 
fj"3"f . Thus in superfluid 4 He, dp c (p)/dp ~ —0.3 for p 
close to pb- In contrast, in a dilute hard sphere Bose 
gas, dp c {p)/dp < 1 §. 

Independently of the numerical results in Refs. |7|,[§, 
it is easy to see that the general structure of the varia- 
tional many-particle wavefunction which LPP and oth- 
ers [p|j9|jr^l6| have used to treat the free surface of liquid 
4 He already builds in the essential physics of a surface re- 
gion which is completely Bose-condensed into one single- 
particle state. In their simplest form, such wavefunctions 
are assumed to have the Feenberg form 



*(ri,...rAr) 



Ae 



[-i$>(**-*<)-4E*( r *>] 



(4) 



including one and two-particle correlations. The function 
t(r) controls the shape of the density profile p(r) associ- 
ated with this wavefunction. It is easy to see that (^) 
can be rewritten in the form 



*(ri, . . .r/v) = A#ja8trow(ri, . . .rjv)* g as(ri, . ..r N ) 



where ^/jastrow alone describes a bulk liquid and 
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*gas(ri 



■ r N j 



(5) 



(6) 



describes a Bose gas of N atoms all occupying the same 
single-particle state given by </>o( r ) = e~2*( r ) (i.e., with 
100% BEC). Choosing u(r) and t(r) appropriately, wave- 
functions such as ([|) allow one the variational freedom 
to describe atoms in the bulk region (where ^gas — 1) 
and in the low density surface region (note that when 
the atoms are far apart, one has ^jastrow — 1). It is 
important to realize that 4>o(r) is not the single-particle 
"natural orbital" in (|J), since the latter is defined relative 
to the complete wavefunction in (||). The interpretation 
of ([|) as described above in (|J) and (JsJ) has not been em- 
phasized in the previous literature on superfluid surfaces 
(see, however, Ref. 0) but we think the present analysis 
clarifies the physics behind the complete BEC found in 
the surface region. 



In the rest of this letter, we discuss how our new picture 
can be used to describe the surface region in a rigorous 
fashion using the field-theoretic analysis of an inhomoge- 
neous Bose-condensed fluid |5U|]. In addition, we point 
out that this has important implications concerning den- 
sity functional approaches used to calculate the surface 
properties of superfluid 4 He [Q,|l7],[l8| . 

The field-theoretical approach to Bose systems is based 
on the presence of the symmetry-breaking order param- 
eter $(r) = (ip(r)) as the anomalous average of the field 
operator ?/>(r) (see Ch. 3 of Ref. ||). The local conden- 
sate density is given by /o c (r) = |$(r)| 2 . The equation 
for the order parameter can be obtained starting from 
the Heisenberg equation for the field operator 



i^(r) = [F',V>(r)] 



(7) 



where H' = H — pN is the grand canonical Hamiltonian. 
By carrying out explicitly the commutator of (Q) and 
taking the statistical average on the equilibrium state of 
the system one obtains the exact equation ]^,[l9| 



h z V 2 
2m 



$(r) 



dr'v(r - r')(^ t (r')V'(r')V'(r)} = , 

(8) 



where v(r — r') is the bare He- He interatomic potential. 
When applied to the homogeneous liquid, (g) provides 
an exact, non-trivial relationship between the chemical 
potential p and the long range behaviour of the non- 
diagonal two-body density matrix of a Bose condensed 
system [fjcfl . 

In the case of the free surface, the value of the chemi- 
cal potential coincides with the bulk 4 He binding energy 
(— 7.15K at T = 0). In this case, (||) can be used to inves- 
tigate the behaviour of the order parameter ^(r) when 
r is in the asymptotic, low density region far from the 
surface. Our reasoning is as follows. The dominant con- 
tribution to the integral in comes from r' in the bulk 
region, where p(r') is large. For this contribution, |r — r'| 
is large and the correlation function in (||) is given by the 
asymptotically exact formula 

(^VMr'M^lir-r'Hoo = $V)^(0>$(r)> =P(r'Mr) 

(9) 

This may be viewed as a generalization of the Penrose- 
Onsager result |2l]] 

(ft(r')i>(r))\ lr _ T , Hoc = $V))$(r)) = $*(r>(r) . 

(10) 

Using ( lit) in (|J), we conclude that for r in the low density 
region, (ph reduces to 

- i" + / d r'v(r - r'Mr')} *(r) - , (11) 
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where the Hartree potential 



v H (r) = / dr'v(r — r')p(r') 



(12) 



describes the field at r due to the bulk region where the 
density is large. For such contributions, where |r — r'| is 
large, only the long-range attractive van der Waals' tail 
of v(r — r') is important. For contributions to (||) when 
|r — r'| is small (i. e. , when r' is in the low surface re- 
gion), one can use (r')i/>(r')?/>(r)) ~ |<E>(r')| 2 <E>(r) since 
we are dealing with a region of complete BEC. Such con- 
tributions are precisely those kept in the standard Gross- 
Pitaevskii theory |2^|| and as usual one needs to include 
additional multiple scattering terms || which screen the 
hard core in v(r — r') in (ph. However the details are not 
important here (see Ref . |p3| ) since the low density at 
r makes these contributions to (^) negligible compared 
to the ones discussed above. Similarly, while the many- 
body screening of the hard core is not explicitly included 



in (y_l|) and (p_2[), contributions to Uij(r) when |r — r'| is 
small are negligible because of the low density at r. 

In summary, we conclude that ( pi] ) allows one to de- 
termine the exact asymptotic behavior of <i>(r) far from 
the surface of superfluid 4 He [^3|. Interpreting ( pi] ) in 
physical terms, we see that the bulk of the liquid pro- 
duces an effective field which acts to stabilize the dilute 
inhomogeneous gas in the surface region. The fact that 
the atoms in the low density region only feel this effective 
field but otherwise are not affected by their interactions 
explains why this region is fully Bose-condensed. Equa- 
tions with the same Hartree-like structure as (|ll]) arise 
in discussions of the binding energy of an impurity atom 
on the free surface of liquid 4 He, based on variational 
wavef unctions [ p4| . 

For a planar surface centered at z — 0, the asymptotic 
behavior of (|TJ) reduces at large z to vh(z) = —pBa/z 3 , 
where a only depends on the parameters of the He-He 
attractive potential p4j| . Using this in ([ll]), one finds 
that the order parameter has the form 



$(z) 



2m | p 



oc e 



-(Az+B/z 2 ) 



(13) 



lA^ 1 and B = m,p B a/2h 2 A 



where A = 

5A 2 . The above discussion then suggests that the ex- 
ponential decay p5[ exhibited by the resulting density 
profile p(z) — p c (z) = |$(z)| 2 is a direct consequence 
of the complete BEC in the low density region, which 
occurs (see Fig. 1) for p(z) < O.lps- This behavior dif- 
fers, for example, from the surface density profile of a 
Fermi liquid, where one finds a faster decay, of the form 
p(z) oc z~ 2 e~ 2Az ||]. Integrating (|l|) to find the total 
number of atoms in the region of essentially 100% BEC, 
it corresponds to an equivalent surface density of about 
10 14 atoms/cm 2 . 

Whatever the density is, the kinetic energy functional 
of a Bose-condensed system always has a contribution 
directly related to $(r) which is given by [fi|,|l9[ 



/ 



dr$*(r)[-^]*(r) = J dv^-\V y/p~Jv)\ 2 

+ i m J drp c (r)v 2 (r) , (14) 

where $(r) = yj p c {v)e lS ^ and mv s (r) = HVS(r) is the 
local superfluid velocity. The expression in (|14| ) is the 
origin of the kinetic energy term in (||) and (p_l|). 

The long-range coherence effects associated with su- 
perfluidity JiiJ are related to the local phase S(r) of the 
complex (two-component) Bose order parameter $(r) de- 
fined above. In a situation like the one considered here, 
where the local condensate density p c (j) is rapidly de- 
creasing normal to the surface, the superfluid flow prop- 
erties may be quite different parallel and perpendicular to 
the surface. This clearly has experimental implications. 

Phenomcnological density functional theories used for 
inhomogeneous superfluid 4 He JT^JT^Q always start from 
an energy functional of p(r) alone, usually of the form 



H[p(r)\ = [ dr^\V^plr)\ 2 + V COII [p(r)} 
J 2m 



(15) 



The first term is the kinetic energy of an inhomogeneous 
ideal Bose gas with a given local density p(r), while 
Korr[p(r)] is a phenomenological expression used to de- 
scribe the effects of interactions. More precisely, the first 
term in (Jlf^) describes the energy of a fictitious liquid of 
N atoms as if t hey all o ccupied the same single-pa rticle 
state <j)Q (r) = y/p(r)/N, corresponding to $(r) = \J p(r). 
A complete density functional theory of Bose-condensed 
liquids would involve functionals of both the local den- 
sity p(r) and the local order parameter 3>(r) 27|. In such 
a theory, the first term o f (1141 ) would naturally arise in 
place of the first term in (Jl 5[) , although we note that the 
two terms coincide in the important low density surface 
region since p c (r) — p(r). 

It would be useful to have finite temperature path in- 
tegral Monte Carlo calculations of p c (r) in the surface re- 
gion of superfluid 4 He [p8| . Physically, it seems clear that 
the surface region will Bose condense at the bulk super- 
fluid transition temperature T\ ~ 2.17K. The chemical 
potential is the same for the liquid and vapor in equilib- 
rium and since it is large and negative (p ~ — 7K from 
to 2K), the 4 He vapor can always be described as an ideal 
classical gas. Fig. [j] is only valid at low temperatures 
(< IK). As the temperature increases, the condensate 
density p c ( r ) will be increasingly depleted ||] and thus 
it will be a smaller component of the total local density 
p(r), in both the bulk and surface regions. However we 
expect that /? c ( r ) w iU still decrease exponentially in the 
surface region, while p(r) will go over smoothly to the 
small but finite density of the 4 He vapor. 

In conclusion, in this letter we have argued that the 
surface region of superfluid 4 He in self-bound phases can 
be described as an inhomogeneous dilute Bose gas which 
is completely Bose-condensed at T = 0. Dramatic evi- 
dence for this "self-bound Bose gas" is given by the nu- 
merical calculations in Ref. 0. We have noted that the 
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asymptotic low density surface region can be described 
exactly by the order parameter $(r) given by the solution 
of a generalized Gross-Pitaevskii equation ( pd| ) taking full 
account of the long range van der Waals tail of the He-He 
potential. Optical excitation of the surface atoms might 
be a way of detecting the completely Bose-condensed na- 
ture of the surface region of superfluid 4 He. We have 
used our new physical ideas to suggest a more micro- 
scopic basis for the kinetic energy term used in density 
functional theories. This inhomogeneous Bose-condensed 
gas is readily available, and should complement the more 
exotic Bose gas systems on which much current research 
has concentrated [ p9[ . 
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